Quantum critical behavior of ultracold atoms in two-dimensional optical lattices 
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As the temperature of a many-body system approaches absolute zero, thermal fluctuations of 
observables cease and quantum fluctuations dominate. Competition between different energies, such 
as kinetic energy, interactions or thermodynamic potentials, can induce a quantum phase transition 
between distinct ground states. Near a continuous quantum phase transition, the many-body system 
is quantum critical, exhibiting scale invariant and universal collective behavior [T] [2]. Quantum 
criticality has been actively pursued in the study of a broad range of novel materials [3HS] , and can 
invoke new insights beyond the Landau-Ginzburg- Wilson paradigm of critical phenomena [7]. It 
remains a challenging task, however, to directly and quantitatively verify predictions of quantum 
criticality in a clean and controlled system. Here we report the observation of quantum critical 
behavior in a two-dimensional Bose gas in optical lattices near the vacuum-to-superfluid quantum 
phase transition. Based on in situ density measurements, we observe universal scaling of the equation 
of state at sufficiently low temperatures, locate the quantum critical point, and determine the critical 
exponents. The universal scaling laws also allow determination of thermodynamic observables. In 
particular, we observe a finite entropy per particle in the critical regime, which only weakly depends 
on the atomic interaction. Our experiment provides a prototypical method to study quantum 
criticality with ultracold atoms, and prepares the essential tools for further study on quantum 
critical dynamics. 

PACS numbers: 05.70.Jk, 64.60.F-, 64.70.Tg, 67.85.Hj 



In the vicinity of a continuous quantum phase transi- 
tion, quantum fluctuations lead to non-classical universal 
behavior of a many-body system pQ. Quantum critical- 
ity not only provides novel routes to new material design 
and discovery pQ, but also raises possible links between 
condensed matter systems and those studied in nuclear 
physics [3 [9] or in cosmology [TJ [TO] . Understanding 
quantum criticality and its general role in strongly cor- 
related systems has hence attracted significant studies 
such as those on heavy- fermion materials [4] [11] , quasi- 
one-dimensional Ising ferromagnets [T2] . quantum anti- 
ferromagnets [13 , ruthenate metals [14], and chromium 
at high pressure [15 . 

Ultracold atoms offer a clean system for a quantitative 
and precise study of quantum phase transitions and crit- 
ical phenomena [TrTH2Q] . As an example, the supernuid- 
to-Mott insulator quantum phase transition is realized by 
loading atomic Bose-Einstein condenstates into optical 
lattices [16]. In recent experiments, universal scaling be- 
havior was observed in interacting Bose gases in three [2T] 
and two dimensions [22] , and in Rydberg gases [23] . Sup- 
pression of the superfluid critical temperature near the 
Mott transition was observed in 3D optical lattices [24] . 
Studying quantum criticality in cold atoms based on 
finite-temperature thermodynamic measurements, how- 
ever, remains challenging and has attracted increasing 
theoretical interest in recent years [25H28] . 

In this letter, we report the observation of quan- 
tum critical behavior of ultracold cesium atoms in a 
two-dimensional (2D) optical lattice across the vacuum- 
to-superfluid transition. This phase transition can be 
viewed as a transition between Mott insulator with zero 
occupation number and superfluid, and can be quan- 



titatively described by the Bose-Hubbard model [29] . 
Our measurement is based on atomic samples near a 
normal-to-superfluid transition, which connects to the 
vacuum-to-superfluid quantum phase transition in the 
zero-temperature limit. At progressively lower temper- 
atures, quantum criticality is revealed in the emergence 
of universal scaling of the equation of state. From the 
equation of state, we extract the quantum critical point 
and the critical exponents, and compare them with the- 
oretical predictions. Furthermore, we observe the break- 
down of quantum criticality at high temperatures and 
estimate an upper thermal energy scale for the quantum 
critical behavior. The derived scaling laws permit a com- 
plete determination of the thermodynamics of the critical 
gas. In particular, we observe a universal, non-zero en- 
tropy per particle in the critical regime, carrying a weak 
dependence on the atomic interactions. 

The quantum phase transition and quantum critical 
regime in this study are illustrated in Fig. [I] The 
zero-temperature vacuum-to-superfluid transition occurs 
when the chemical potential \i approaches the quantum 
critical point /io- Sufficiently close to the quantum criti- 
cal point, the critical temperature T c for the normal-to- 
superfluid transition is expected to decrease according to 
the following scaling [29] : 

— = c \-^) ' (1) 

where &b is the Boltzmann constant, t is the tunneling 
energy, z is the dynamical critical exponent, v is the cor- 
relation length exponent, and c is a constant. In the 
quantum critical regime (shaded area in Fig. [I]), the tem- 
perature T provides the sole energy scale, and all ther- 
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FIG. 1: The vacuum-to-superfluid quantum phase 
transition in 2D optical lattices. At zero temperature, 
a quantum phase transition from vacuum (horizontal thick 
blue line) to superfluid occurs when the chemical potential \i 
reaches the critical value fio. Sufficiently close to the transi- 
tion point /io, quantum criticality prevails (red shaded area), 
and the normal-to-superfluid transition temperature T c (mea- 
surements shown as empty circles; see Methods) is expected 
to vanish as T c ~ (/i — \±§) zv \ the blue line is a guide to the eye. 
From the prediction zv — 1 [26] |27l [29] , the linearly extrap- 
olated critical chemical potential is fio — — 3.6(6)t, consistent 
with the theoretical value — 4t [18] . Here both the thermal 
energy scale k&T and chemical potential \i are normalized by 
the tunneling t. 



modynamic observables are expected to scale universally 
with T [29]. Thus the equation of state is predicted to 
obey the following scaling [25H27] 

N = F(fi), (2) 
in which F(x) is a universal function, and 

N — N H — m - 

are the scaled occupation number and scaled chemical 
potential, respectively. Here N is the occupation number, 
D is the dimensionality, and N r is the non-universal part 
of the occupation number. For the vacuum-to-superfluid 
transition in the two-dimensional Bose-Hubbard model, 
we have N r = and D = 2, and the predicted critical 
exponents are z = 2 and v — 1/2, characteristics of the 
dilute Bose gas universality class [2j [26j [29] . 

Our experiment is based on 2D atomic gases of cesium- 
133 in 2D square optical lattices [30] [31]. The 2D trap 
geometry is provided by the weak horizontal (r-) confine- 
ment and strong vertical (z-) confinement, with envelope 
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FIG. 2: Evidence of a quantum critical regime, a, 

Scaled occupation number TV = Nt/k^T as a function of the 
scaled chemical potential jl — (fi — /io)/feT, measured at 
seven temperatures: T =5.8 nK (black circles), 6.7 nK (red 
triangles), 11 nK (green triangles), 13 nK (blue diamonds), 
15 nK (magenta squares), 24 nK (orange circles), and 31 nK 
(dark yellow triangles), with the blue solid line showing the 
average curve for the lowest four temperatures. Inset shows 
the low-temperature data in the range of T — 5.8 ~ 15 nK, 
and the critical chemical potential jio is identified as the cross- 
ing point; see text. The result, /j,o — — 4.5(6)£ agrees with the 
prediction fio = —At [18]. b, Determination of the dynam- 
ical critical exponent z and the correlation length exponent 
v. The color represents the sum of normalized mean square 
deviations of N for data in the low temperature range of 
T — 5.8 ~ 15 nK, and indicates how well the scaled equa- 
tion of state can collapse into one single curve (smaller devi- 
ations suggest better fits); see text. Best fit corresponds to 
z = 2.0(3) and v — 0.54(5). The uncertainties are determined 
from the contour on the z-v plane where the deviation is 50% 
above the minimum. The predicted values are z — 2 and 
i/ = 1/2. 



trap frequencies f r = 9.6 Hz and f z = 1940 Hz, respec- 
tively. Details of the system and the sample preparation 
are described in Ref. [31]. Typically 4,000 to 20,000 
Bose-condensed atoms are loaded into the lattice. The 
lattice constant is d = A/2 = 0.532 /im and the depth 
is Vl = 6.8 £?r, where = ^ x 63.6 nK is the re- 
coil energy, A = 1064 nm is the lattice laser wavelength, 
and h is the Planck constant. In the lattice, the tunnel- 
ing energy is t = &b x 2.7 nK, the on-site interaction is 
U = &b x 15 nK, and the scattering length is a = 15.9 nm. 
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FIG. 3: Finite-temperature effect on quantum criti- 
cal scaling. Scaled occupation number iV c = N c t/k B T at 
the critical chemical potential /i = fio as a function of the 
normalized temperature h&T/t. The blue dashed line is an 
empirical fit, giving a temperature scale T* « 8£/&b. For 
T < T*, 7V C ~ 0.097 is independent of the temperature; for 
T > T*, iV c deviates from the low-temperature value. 



The sample temperature is controlled in the range of T 
= 5.8-31 nK. 

We determine the equation of state n(/i, T) of the 
sample from the measured in situ density distribution 
n(x,y) [22, 30]. The chemical potential ja(x^y) and the 
temperature T are obtained by fitting the low-density 
tail of the sample where the atoms are normal. The fit 
is based on a mean-field model which accounts for inter- 
action [18j|32]; see Methods. 

We locate the quantum critical point by noting that at 
the critical chemical potential \i = /io, the scaled occupa- 
tion number N = Nt/k^T = nd 2 t/k&T is temperature- 
independent, as indicated by Eq.JS] We plot N as a func- 
tion of /i in the low temperature range of 5.8—15 nK, and 
indeed observe a crossing point at /io = — 4.5(6)£, see in- 
set of Fig. [2^i. We identify this point as the critical point 
for the vacuum-to-superfluid transition, and our result 
agrees with the prediction —At [18] . 

To test the critical scaling law, we compare the equa- 
tion of state at different temperatures. Based on the ex- 
pected exponents z = 2 and v — 1/2, we plot the scaled 
occupation number TV as a function of the scaled chemical 
potential fl = (/i — ii$)/k&T] see Fig. [2^. Below 15 nK, all 
the measurements collapse into a single curve, which con- 
firms the emergence of the quantum critical scaling law 
(Eq. [2| at low temperatures. Deviations become obvious 
at higher temperatures. 

We determine the critical exponents z and v by com- 
paring the equation of state at low temperatures. Taking 
various values of z and v in the range of < z < 4 and 
< v < 1, we compute the corresponding scaled oc- 
cupation numbers N and scaled chemical potentials jl 
based on Eq. [3] We then evaluate how well the scaled 
equation of state in the range of T = 5.8 — 15 nK can 
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FIG. 4: Entropy per particle in the critical regime, a, 

Entropy per particle S/N as a function of the scaled chemical 
potential jl, measured in the temperature range of 5.8 ~ 15nK 
(same symbol and color scheme as in Fig. [2^i). b, Critical 
entropy per particle S c /N as a function of the effective in- 
teraction strengths g: measurements for Bose gases with 2D 
optical lattices (black circle) and without lattice (black tri- 
angles, extracted from data in Ref. [22]), mean- field calcula- 
tions (blue line), and a power- law fit to the measurements, 
S c /Nk B = 1.6(l)#°- 18(2) (red line). 



collapse to a single curve by computing the normalized 
mean-square deviations of TV (Methods). The result is 
shown in Fig. [2}}, from which we find the exponents to 
be z — 2.0(3) and v = 0.54(5). Our result is consistent 
with the predictions of the dilute Bose gas universality 
class: z = 2 and v = 1/2 [2}[26| 29 , and shows the power 
of our method in identifying universality classes. In the 
following analyses, we adopt z = 2 and v = 1/2. 

Our measurements at different temperatures allow us 
to investigate the breakdown of quantum criticality at 
high temperatures. To quantify the deviations, we fo- 
cus on the temperature dependence of the scaled oc- 
cupation number N at the critical chemical potential 
\i = /io, as shown in Fig. [3] Deviations from the low- 
temperature value are clear when the temperature ex- 
ceeds T* = 22 nK^ St/ks- From this, we conclude that 
at \i = /io, the upper bound of thermal energy for the 
quantum critical behavior in our system is k^T* w St. 
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Our result is in fair agreement with the prediction of 6t 
based on quantum Monte Carlo calculations [27] . 

From the equation of state, one can derive other ther- 
modynamic quantities in the critical regime. Here we de- 
rive the entropy per particle S/N based on measurements 
in the temperature range ofT = 5.8 ~ 15 nK, using a 
procedure similar to that in Ref. [33]. The measured 
entropy per particle only depends on the scaled chem- 
ical potential fl and monotonically decreases (Fig. |4|i), 
indicating a positive specific heat. Near the critical 
point jl = 0, the entropy per particle has an approxi- 
mate linear dependence on the scaled chemical potential: 
S/Nk B = a -bp,, with a = 1.8(1), b = 1.1(1). From 
this linear dependence, we derive an empirical thermody- 
namic relation analogous to the ideal gas law (Methods) : 

P = Cn x (k B T) y , (4) 

where P is the pressure of the 2D gas, x = = 0.95(5), 
y = ^ = 1.05(5), C = 0.8(2)(td 2 ) w is a constant, and 

w = TT5 = -0-05(5). 

Finally, we observe a weak dependence of the criti- 
cal entropy per particle on the atomic interaction. Not- 
ing that a weakly-interacting 2D Bose gas follows sim- 
ilar scaling laws near fi = [22 due to the same un- 
derlying dilute Bose gas universality class [2j [34], we 
apply similar analysis and extract the critical entropy 
per particle S c /N at four interaction strengths g « 
0.05,0.13,0.19,0.26, as shown together with the lattice 
data (g ~ 2.4) in Fig. |4Jd. We observe a slow growing of 
S c /N with g, and compare the measurements with mean- 
field calculations. The measured S/N is systematically 
lower than the mean-field predictions, potentially due to 



quantum critical physics. The weak dependence on the 
interaction strength can be captured by a power-law fit 
to the data as S c /Nk B = 1.6(1)#°- 18 ( 2 ). 

In summary, based on in situ density measurements 
of Bose gases in 2D optical lattices, we confirm the 
quantum criticality near the vacuum-to-superfluid 
quantum phase transition. We show the suppression of 
the superfluid critical temperature, observe the universal 
scaling of equation of state, and extract the quantum 
critical point and critical exponents. In addition, we find 
that the entropy per particle is temperature-independent 
in the critical regime and has a weak dependence on 
the atomic interaction. Our experimental methods hold 
promise for identifying other quantum phase transitions, 
and prepare the tools for investigating quantum critical 
dynamics. 



Methods Summary 

Preparation and detection of cesium 2D Bose gases 
in optical lattices are similar to those described in 
Refs. [22j[3T]. We obtain the normal-to-superfluid tran- 
sition point by collapsing the scaled equation of state 
at different temperatures [22]. We extract the quantum 
critical points and critical exponents by minimizing the 
normalized mean square deviations of scaled occupation 
number. Based on the Gibbs-Duham equation, we derive 
thermodynamic quantities, in particular, the entropy per 
particle in the critical regime (details in Online Meth- 
ods). 
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Methods 

Preparation and characterization of cesium 2D 
Bose gases in optical lattices 

Preparation and detection of 2D gases in optical lattices 
The experimental procedure is similar to those described 
in Refs. [221 El] . We adjust the atomic temperature by 
applying magnetic field pulses near a magnetic Feshbach 
resonance [35] to excite the atoms [22 . After the pulse, 
we tune the scattering length to a = 15.9 nm, wait for 
200 ms, and ramp on the optical lattice to 6.8 in 270 
ms. These parameters are chosen to allow the sample to 
reach thermal equilibrium after the ramp [51] . The final 
lattice depth is sufficiently deep to validate single band 
Bose-Hubbard description. After preparing the sample, 
we perform in situ absorption imaging using a strong 
resonant laser beam [22] [33] . The atomic density is in- 
dependently calibrated in similar methods as in Ref. [22] . 

Determination of the peak chemical potential and the 
temperature We determine the peak chemical potential 
/i m (chemical potential at the trap center r = 0) and the 
temperature T of a 2D Bose gas in 2D optical lattices by 
fitting the low-density tail of the azimuthally averaged 
density profile using the following formula [T5] : 

CO 

n(r) = d~ 2 ^[I (2ipt)] 2 e l ^- 2U ^ nd2 - v ^\ (5) 
i=i 

This formula is based on local density approximation and 
a mean-field model that takes interaction into account. 
Here n(r) is the 2D atomic density at radius r from 
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the cloud center, d = 0.532 /im is the lattice spacing, 
Io(%) = fin ^e xcos9 is the zeroth-order Bessel function 
with purely imaginary argument, f3 = l/k B T, k B is the 
Boltzmann constant, t is the tunneling, V(r) is the en- 
velope trapping potential, and C7 e ff is the effective inter- 
action. Here the calculation of U e R involves the Bose- 
Hubbard on-site interaction parameter U and terms for 
a modified two-particle propagator [32] : 



Z7 eff = 



U 



where 

n = 



2t LL 



d/c x d/c v 



(6) 



fc B T 
t 



2 [2 — cos(k x d) — cos(kyd)} ' 



and the integration ranges of fc x and k y are both from 
—ir/d to n/d, which covers the first Brillouin zone of 
the 2D square lattice. We test this formula on quantum 
Monte Carlo (QMC) data [27]. Within our experimental 
temperature range, the fitted T agrees with QMC value 
within 3%, and the fitted /i m agrees with QMC value 
within 0.6t. 

Determination of the critical parameters 

Normal-to- super fluid transition point /i c We use a pro- 
cedure similar to that in Ref. [22 . At a reference tem- 
perature T re f, we obtain the critical chemical potential 
/i c ,ref for the normal-to-superfluid transition by fitting 
the crossover of the compressibility k = ^ (as a function 
of density n) near the transition region. The fit is based 
on an empirical formula k = kn — ^k 2 {n — n c ) 2 + w 2 + 
^k 2 n 2 + w 2 , where the slope k, critical density n c = 
n(/z re f , T re f ), and the width w of the transition region are 
fitting parameters. At a different temperature T, we ob- 
tain the transition point by comparing the equation of 
state at the two temperatures T and T re f and collapsing 
the scaled equation of state near the transition points 
into a single universal curve [22] : 



k B T k B T 
The resulting critical points are shown in Fig. [T] 



(7) 



Quantum critical point /aq For a given chemical potential 
/i, we calculate the mean square deviation of the scaled 
occupation numbers N = Nt/k B T at M different tem- 
peratures, and normalize this deviation by the squared 
average value, as shown in the following formula: 



1 



{M-l)Nl{ti)U 



M 2 

J2[Ni(»)-NM] (8) 



where 7V av (/i) = YliL\Ni{lj) /M . The quantum critical 
point /i = fiQ is determined by finding the minimum of 
AAT(/i). 



Critical exponents z and v For given trial values of z and 
v, we calculate the scaled occupation number N and 
the scaled chemical potential ft according to Eq. [3] and 
define the normalized mean square deviation, AiV(/i), at 
a certain jl in a form similar to Eq. [8] The exponents z 
and v are determined by minimizing the average AN(ji) 
in the range of —1.5 k B T < /i — fiQ < 1.5 /cbT. 

Universal thermodynamics in the quantum crit- 
ical regime Based on the Gibbs-Duham equation [36] . 
we derive the pressure P(/i, T) from the in situ density 
measurements: 



P(fi,T) 



n(fj,' ,T)d/j? 



(9) 



The entropy density s(fi, T) is related to the pressure via 
differentiation with respect to temperature: 



zap 

V dT 



(10) 



In the quantum critical regime near the vacuum-to- 
superfluid transition, the scale invariance of the equation 
of state suggests the following scaling laws for pressure 
P and entropy density s: 



P(JJL,T) 



s(fi,T) 
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K P (ji) , 



K s (A) , 



(11) 



(12) 



where Kp and K s are generic functions of jl. Combining 
Eq.[2]and Eq. 12, we obtain the entropy per particle S/N 
in unit of k B 



Nk ¥ 



k B T 
t 



-i+- 



(13) 



where W is a universal function. 

For z = 2, v = 1/2, S/Nk B = W{p) = 2K P /K P ' - 
jl is a temperature-independent function of the scaled 
chemical potential jl = (ji — jio) / k B T . In particular, near 
the critical point jl = 0, S/Nk B varies approximately as 
a linear function of jl: W{jl) = a — bjl. Using W{jl) = 
2Kp/Kp' — jl, we can solve the pressure P and express 
it in terms of density n and temperature T: 



where x 



2 

1+6' 



Cn x {k B T)y, (14) 
and the proportionality con- 



A similar technique was applied to obtain the entropy 
per particle for a bulk 2D gas of rubidium-87 atoms [33] . 



stant C = (f) T +^(- 
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Effective interaction strength of a 2D gas We define the 
dimensionless effective interaction strength g for our Bose 
gas in optical lattices: 



9 = U e f£- 



m*d 2 



and for that without lattices 1221: 



(15) 



oscillator length. 

Mean-field calculation on the entropy per particle We cal- 
culate the entropy per particle S/N based on Eq. [5] At 
low temperatures T, the Bessel function takes its asymp- 
totic form Iq(x) ~ e x '/y/2irx when x = 2lt/{k^T) is large, 
and Eq. [5] reduces to 



(16) 



where the effective 
using Eq. § h ' 



interaction is calculated 

is the reduced Planck constant, 
vrC = H 2 /E"(k)\k = Q is the single-particle effective 
mass in a 2D optical lattice and can be calculated 
from the ground-band dispersion relation E(k), a is 
the scattering length which is tunable via a magnetic 
Feshbach resonance [35], and L is the vertical harmonic 



In 



exp [fi —F(fi) 



One can calculate F{p) by solving this equation self- 
consistently, and then derive Kp(jl) and S/Nkft = 
2Kp I 'Kp' —jl from F(p). In this calculation, the effective 
mass takes the value m* = (under the tight-binding 
approximation); the effective interaction strength is thus 
given by g = -^f- based on Eq. 15 



